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Abstract 

We investigate deep inelastic scattering from 3 He and 3 H within a conven- 
tional convolution treatment of binding and Fermi motion effects. Using re- 
alistic Faddeev wave functions together with a nucleon spectral function, we 
demonstrate that the free neutron structure function can be extracted in deep- 
inelastic scattering from A = 3 mirror nuclei, with nuclear effects canceling 
to within 2% for x < 0.85. 
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One of the most fundamental properties of the nucleon is the structure of its valence quark 
distributions. Unlike the sea, which is generated via both perturbative and non-perturbative 
mechanisms, the valence quark structure reflects entirely large distance dynamics in the 
nucleon, which cannot be described within perturbative quantum chromodynamics. 

Experimentally, most of the recent studies of nucleon structure have emphasized the 
small-x region populated mainly by sea quarks (x being the fraction of momentum of the 
nucleon carried by the quark) , while the valence quark structure has for some time now been 
thought to be understood. This is to some extent true, albeit with one major exception 
- the so-called deep valence region, at very large x, x > 0.7. Recently it has become 
more widely appreciated that knowledge of quark distributions at large x is essential for 
a number of reasons. Not least of these is the necessity of understanding backgrounds in 
collider experiments, such as in searches for new physics beyond the standard model [|TJ. 
Furthermore, the behavior of the d/u quark distribution ratio in the limit x — > 1 is a critical 
test of the mechanism of spin-flavor symmetry breaking in the nucleon, and of the onset of 
perturbative behavior in large- £ structure functions. 

The most widely used source of information about the valence quark distributions in the 
nucleon has been the proton structure function, Ff, which at large x measures a charge- 
squared weighted combination of the valence u and d distributions. Because the u quark 
is weighted by a factor 4:1 compared with the d, the structure function most directly 
constrains the u quark distribution. 

To determine the individual isospin distributions separately requires a second linear 
combination of u and d, which traditionally is obtained from the neutron structure function, 
and which could in principle constrain the d quark distribution as well as the u. However, 
the absence of free neutron targets means that usually the deuteron is used as an effective 
neutron, with the neutron structure function approximated by F£ ~ F 2 rf - F$. While this 
approximation is valid at moderate x, it breaks down dramatically for x ^ 0.4 due to Fermi 
motion and nuclear binding effects in the deuteron . 

The problem of extracting neutron structure functions from nuclear data is rather old 
0, although recently the discussion has been revived with the realization || that F£, 
extracted from F^ by taking into account Fermi motion and binding (off-shell) effects, could 
be significantly larger than that extracted in earlier analyses in which only Fermi motion 
corrections were applied. In particular, omitting nuclear binding corrections can introduce 
errors of up to 50% in F% jF% already at x ~ 0.75. Such a difference is of the same order 
of magnitude as the variation of the behavior of the F^jF^ ratio predicted in the x — > 1 
limit, which ranges from 1/4 in non-perturbative models where the d quark is suppressed 
relative to u H, to 3/7 in perturbative QCD- inspired models which emphasize helicity 
aligned configurations of the quark and nucleon 0. 

Although one can make a strong argument that a proper treatment of nuclear corrections 
in the deuteron should account for both Fermi motion as well as binding effects, the question 
can ultimately be settled only by experiment. In this paper we suggest how this can be 
achieved by using a novel method which maximally exploits the mirror symmetry of A = 3 
nuclei. Regardless of the absolute value of the nuclear EMC effects in 3 He or 3 H, the 
differences between these will be small — on the scale of charge symmetry breaking in 
the nucleus — which allows a relatively clean determination of F% over a large range of x 
essentially free of nuclear contamination. 
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The argument is actually rather simple. In the absence of the Coulomb interaction and 
in an isospin symmetric world the properties of a proton (neutron) bound in a 3 He nucleus 
would be identical to that of a neutron (proton) bound in 3 H. If, in addition, the proton 
and neutron distributions in 3 He (and in 3 H) were identical, the neutron structure function 
could be extracted with no nuclear corrections, regardless of the size of the EMC effect in 
3 He or 3 H separately. 

In practice, 3 He and 3 H are of course not perfect mirror nuclei — their binding energies 
for instance differ by some 10% — and the p and n distributions are not quite identical. 
However, the A = 3 system has been studied for many years, and modern realistic A = 3 
wave functions are known to rather good accuracy. In a self-consistent framework one can 
use the same NN interaction which describes the two-nucleon system (e.g. NN scattering, 
deuteron form factors, quasi-elastic ed scattering) to provide the basic input interaction into 
the three-nucleon calculation. Therefore the wave functions can be tested against a large 
array of observables which put rather strong constraints on the models. 

Defining the EMC-type ratios for the F 2 structure functions of 3 He and 3 H (weighted by 
corresponding isospin factors) by: 

Ft 



* (3He) = WW? • (la) 

F$ + 2F? 



2FI + F% 

^( 3 H) = J\™ > ( lb ) 



one can write the ratio of these as: 

^ i?( 3 He) 

Inverting this expression directly yields the ratio of the free neutron to proton structure 
functions: 



F 2 n _ 2K - F 2 Re /F 2 R 
F! ~ 2F 2 3He /F 2 3H - TZ 



(3) 



We stress that F 2 / F 2 extracted via Eq. @ does not depend on the size of the EMC effect 
in 3 He or 3 H, but rather on the ratio of the EMC effects in 3 He and 3 H. If the neutron and 
proton distributions in the A = 3 nuclei are not dramatically different, one might expect 
1Z « 1. To test whether this is indeed the case requires an explicit calculation of the EMC 
effect in the A = 3 system. 

The conventional approach employed in calculating nuclear structure functions in the 
valence quark region, x > 0.3, is the impulse approximation, in which the virtual photon 
scatters incoherently from individual nucleons in the nucleus ||. The nuclear cross section 
is determined by factorizing the 7*-nucleus interaction into 7*-nucleon and nucleon-nucleus 
amplitudes. In the absence of relativistic and nucleon off-shell corrections [[T^JXTj] , the struc- 
ture function of a nucleus, F 2 A , can then be calculated by folding the nucleon structure 
function, F 2 , with a nucleon momentum distribution in the nucleus, Jn/a- 

F 2 A (x) = Jdy f(y) F 2 N (x/y) = f(x) ® F 2 N (x) , (4) 
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where y is the fraction of the 'plus '-component of the nuclear momentum carried by the 
interacting nucleon, and the Q 2 dependence in the structure functions is implicit. The 
convolution expression in Eq.(|4]) is correct in the limit of large Q 2 ; at finite Q 2 there are 
additional contributions to F 2 A from the nucleon structure function, although these are 
suppressed by powers of M 2 /Q 2 , where M is the nucleon mass. Corrections to the impulse 
approximation appear in the guise of final state interactions, multiple rescattering (nuclear 
shadowing), NN correlations and 6-quark clusters, however, these are generally confined to 
either the small-x fli~2|| , or very large-x (x ^ 0.9) |13| regions. 

The distribution f(y) of nucleons in the nucleus is related to the nucleon spectral function 
Sip) by O: 




P0+Pz 

M 



Sip) , (5) 



where p is the four-momentum of the bound nucleon, and is normalized such that / dy f(y) = 
1. The spectral function includes final state interactions between the two spectator nucleons, 
either as a deuteron or in the continuum. In the latter case, this includes an integration 
over interacting final state NN wave functions, as distinct from a calculation in terms of 
a simple momentum distribution (see Ref. |14| for a definition of the three-body spectral 
function). For an A = 3 nucleus the spectral function is evaluated from the three-body 
nuclear wave function, calculated by solving the homogeneous Faddeev equation with a 
given two-body interaction. Details of the computation of the wave functions can be found 
in Ref. ||15|| . To examine the model dependence of the distribution function we use several 



different potentials, namely the "EST" (Ernst-Shakin-Thaler) separable approximation to 



the Paris potential |T6[ (referred to as "PEST"), the unitary pole approximation [17] to the 



Reid Soft Core (RSC) potential, and the Yamaguchi potential |TJ| with 7% mixing between 
3 Si and 3 D 1 waves. 

In terms of the proton and neutron momentum distributions, the F2 structure function 
for 3 He is given by: 

F 2 He = 2 / p/3He <g> Ff + f n/me ® F 2 n . (6a) 

Similarly for 3 H, the structure function is evaluated from the proton and neutron momentum 
distributions in 3 H: 

i£ H = / P /»H ® F$ + 2 / n/ 3 H ® F 2 " . (6b) 

Because isospin symmetry breaking effects in nuclei are quite small, one can to a good 
approximation relate the proton and neutron distributions in 3 He to those in 3 H: 

fn/m ~ /p/3He , (7a) 
/p/3H ~ fnpKe , (7b) 

although in practice we consider both the isospin symmetric and isospin symmetry breaking 
cases explicitly. Note that even in the isospin symmetric case the proton and neutron 
distributions in 3 He will be different because while the neutron in 3 He is accompanied by 
a spectator pp, the spectator system of the proton is either an uncorrelated pn pair or a 
recoiling deuteron. 
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The ratio TZ of EMC ratios for He and H is shown in Fig. 1 for the various nuclear 



model wave functions (PEST, RSC and Yamaguchi), using the CTEQ parameterization [19 
of parton distributions at Q 2 = 10 GeV 2 for . The EMC effects are seen to largely cancel 
over a large range of x, out to x ~ 0.85 — 0.9, with the deviation from a 'central value' 
TZ ~ 1.01 within ±1%. Furthermore, the dependence on the nuclear wave function is very 
weak. In practice, the exact shape of TZ will not be important for the purposes of extracting 
F 2 /F 2 from the F 2 He /F 2 H ratio; rather, it is essential that, as we find, the model dependent 
deviation of TZ from the central value should be small. 

The dependence of TZ on the input nucleon structure function parameterization is il- 
lustrated in Fig. 2, where several representative curves at Q 2 = 10 GeV 2 are given: apart 
from the standard CTEQ fit (solid), the results for the GRV [^(| (dot-dashed), Donnachie- 
Landshoff (DL) |2T| (dashed), and BBS |22| (dotted) parameterizations are also shown (the 
latter at Q 2 = 4 GeV 2 ). For x < 0.6 there is little dependence (< 0.5%) in the ratio on 
the structure function input. For 0.6 < x < 0.85 the dependence is greater, but still with 
< ±1% deviation away from the central value TZ ~ 1.01. The spread in this region is due 
mainly to the poor knowledge of the neutron structure function at large x. Beyond x « 0.85 
there are few data in the deep-inelastic region on either the neutron or proton structure 
functions, so here both the d and u quark distributions are poorly determined. 

A standard assumption in most global fits of parton distributions is that d/u — > as 
x — > 1. This assumption has recently been questioned on theoretical and phenomenological 
grounds P, p3p4 |. The BBS parameterization | 22]| , on the other hand, incorporates con- 
straints from perturbative QCD, and forces d/u — *> 0.2 as x — > 1 ||. The effect of the 
different large-x behavior of the d quark is apparent only for x ^ 0.85, where it gives a dif- 
ference of ~ 1-2% in TZ compared with the fits in which d/u — > 0. One can also modify the 
standard CTEQ fit, for example, by applying a correction factor K3] to enforce d/u — > 0.2, 



however, this also produces differences in TZ which are ^ 2% for x < 0.9. 

Despite the seemingly strong dependence on the nucleon structure function input at 
very large x, this dependence is actually artificial. In practice, once the ratio F 2 He /F 2 u is 
measured, one can employ an iterative procedure to eliminate this dependence altogether. 
Namely, after extracting F 2 /F 2 from the data using some calculated TZ, the extracted F 2 
can then be used to compute a new TZ, which is then used to extract a new and better value 
of F 2 /F 2 . This procedure is iterated until convergence is achieved and a self-consistent 
solution for the extracted F 2 / F 2 and TZ is obtained (see also Ref. fl25|). 

All of the structure functions discussed thus far have been calculated assuming leading 
twist dominance at Q 2 = 10 GeV 2 . To test the sensitivity of the ratio to possible effects 
beyond leading twist, we have calculated TZ using the fit to the total F 2 structure function 
from Donnachie and Landshoff [21], which has an explicit higher twist (oc 1/Q 2 ) component 
in addition to the leading twist. The result is indicated by the upper dot-dashed curve 
DL(HT) in Fig. 2. The difference between the leading twist only and leading + higher twist 
curves is negligible for x < 0.8, increasing to ~ 1.5% at x ~ 0.85, where higher twist effects 
are known to be more important. The size of the higher twist corrections can be determined 
by taking measurements at several values of Q 2 and observing any 1/Q 2 dependence of the 
structure function. In particular, since the Q 2 dependence of F 2 has been measured in a 
number of earlier experiments [26], the Q 2 dependence of the extracted F 2 /F 2 can be used 
to separate the leading twist from the non-leading twist components of F 2 . 



5 



We conclude therefore that the effect on TZ from the present lack of knowledge of the 
nucleon structure function is ^ 2% for x ^ 0.85. However, this uncertainty can in principle 
be eliminated altogether via an iteration procedure, so that the only model dependence of 
TZ will be from the nuclear interaction in the A = 3 nucleus. 

The ratios in Fig. 1 were calculated using three-nucleon wave functions neglecting the 
Coulomb interaction and working in an isospin basis (we also omit possible three-body forces 
since these are expected to have a negligible effect on TZ) . To estimate the effect of neglecting 
the Coulomb interaction in 3 He and at the same time correct the long range part of the 
three-body wave function due to the change in the binding energy, we have modified the 1 5'o 
potential in 3 He and 3 H to reproduce their respective experimental binding energies. In this 
way the 3 Si — 3 D\ interaction responsible for the formation of the deuteron is unchanged. 
This approximation spreads the effect of the Coulomb interaction over both the pp and np 
interaction in the 1 So channel. To that extent, it shifts some of the Coulomb effects in the 
neutron distribution in 3 He to the proton distribution. However, this simple modification to 
the 1 S interaction will allow us to study explicitly the possible effects associated with the 
differences in the binding energies of 3 He and 3 H. 

The ratio TZ calculated with the PEST wave function modified according to this prescrip- 
tion is shown in Fig. 3, labeled PEST(E) (dashed curve). (The CTEQ parameterization of 
the nucleon structure function at Q 2 = 10 GeV 2 is used.) The result of this modification 
is a shift of approximately 0.5-1% shift in TZ, with the net effect still being a ratio which 
deviates by < 2% from unity. 

Also shown in Fig. 3 is the prediction of the nuclear density model, extrapolated from 
heavy nuclei to A = 3 |[27|| . The nuclear density model, which has proven successful for 



studying the A-dependence of the EMC effect for heavy nuclei, stems from the empirical 
observation that for heavy nuclei the deviation from unity in the EMC ratio is assumed to 
scale with nuclear density: 

R(A 2 ) - 1 p(A 2 ) ' 1 ) 

where p(A) is the mean nuclear density. From the empirical A = 3 charge radii one finds 
that p( 3 H)/p( 3 He) » 140%, so that the EMC effect in 3 H is predicted to be ~ 40% bigger 
than in 3 He. However, assuming that i?( 3 He) can be extrapolated from the measured EMC 
ratios for heavy nuclei such as 56 Fe, one still finds that ratio \1Z — 1| < 2% for all x < 0.85. 
Although there are questions about the meaning of nuclear density for a few-body system 
28[, it is reassuring to see that practically the entire range of models of the nuclear EMC 



effect predict that TZ is within 1-2% of unity for all x % 0.85. 

The ideal place to carry out a high-x deep-inelastic scattering (DIS) experiment on 3 He 
and 3 H |29|j30[1 is Jefferson Lab ( JLab) with its proposed energy upgrade to 12 GeV. Since the 
ratio of longitudinal to transverse photoabsorption cross sections R — glI^t is the same for 
3 He and 3 H, measurements of the 3 He and 3 H DIS cross sections under identical conditions 
can provide a direct measurement of the ratio of the F 2 structure functions of the two nuclei: 
cr( 3 H)/cr( 3 He) = F^/F^. The key issue for the experiment will be the availability of a 
high density 3 H tritium target similar to those used in the past to measure the elastic form 
factors of 3 H at Saclay |JT] and MIT-Bates ||32|| . The high intensity of the JLab beam and 



the large acceptance of existing or proposed JLab spectrometers will facilitate high statistics 
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DIS cross section measurements (< ±0.25%) over a large x range (0.10 < x < 0.83) and 
valuable systematics checks in a data taking period of just a few weeks. 

The measured F 2 H /F 2 He ratio is expected to be dominated by experimental uncertainties 
that do not cancel in the DIS cross section ratio of 3 H to 3 He, and the theoretical uncertainty 
in the calculation of TZ. Assuming that the target densities can be known to the ~ 0.5% level 
and that the relative difference in the 3 H and 3 He radiative corrections would be ~ 0.5%, 
the total experimental error in the DIS cross section ratio of 3 H to 3 He should be < 1.0% 
(similar to the error of past DIS measurements of the proton to deuteron cross section ratio 
J3Tj|). Such an error is comparable to the present theoretical uncertainty in the calculation 
of the ratio TZ. 

Figure 4 shows the presently available data on F2/F2, adjusted for the JLab 12 GeV 
kinematics, as extracted from the SLAC deep-inelastic a(p) and a(d) cross sections using a 
Fermi-smearing model with the Paris nucleon-nucleon potential || . To indicate the quality 
of the proposed F£ /F^ ratio determination from the cr( 3 H)/o"( 3 He) measurement, we plot in 
Fig. 4 the ± one standard deviation projected error band for the x range accessible with a 12 
GeV upgraded JLab beam. The band includes both projected experimental and theoretical 
uncertainties. The central values of the band represent F^/Fr? determined using the density 
model J27I] for the nuclear EMC effect and data on EMC ratios for heavy nuclei from the 
SLAC experiment E139 fl3"4" |. It is evident, therefore, that the proposed measurement will 
be able to unambiguously distinguish between the two different methods of extracting the 
F2/F2 ratio from proton and deuterium DIS measurements, and determine its value for 
large x with an excellent precision in an (almost) independent model way. 

As well as offering a relatively clean way to extract F2/F2 , DIS from the 3 He/ 3 H system 
can also determine the absolute size of the EMC effect in A < 3 nuclei. With F 2 n determined 
from the combined F 2 He /.F 2 H and F 2 structure functions, the size of the EMC effect in the 
deuteron (namely, F^jiF^ + -^2")) can be deduced from the measured i^/Ff ratio. This 
would settle a question which has remained controversial since the early 1970s. Furthermore, 
data on the absolute values of F 2 He and F 2 H will also allow the absolute value of the EMC 
effect in A = 3 nuclei to be determined. To date the only data on -F 2 He in existence are 
those from the HERMES experiment which measured the ratio a( 3 He)/(a(d) + er(p)), 
although the focus there was the region of small x and Q 2 . 

In summary, we have demonstrated the effectiveness of using A = 3 mirror nuclei to 
extract the ratio of the neutron to proton structure functions, F^/F^, essentially free of 
nuclear effects for all x < 0.85. A successful program of DIS measurements of A = 3 cross 
sections at an energy-upgraded Jefferson Lab would not only settle a "text-book" issue which 
has eluded a definitive resolution for nearly 30 years, but would also allow the completion of 
the empirical study of nuclear effects in deep-inelastic scattering over the full range of mass 
numbers. 
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FIG. 1. Ratio of nuclear EMC ratios for 3 He and 3 H for various nuclear models: PEST (solid), 
Reid Soft Core (dashed), Yamaguchi (dot-dashed). 
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FIG. 2. Ratio of nuclear EMC ratios for 3 He and 3 H with the PEST wave functions, using 
various nucleon structure function parameterizations: CTEQ (solid), GRV (dot-dashed), BBS 
(dotted), and Donnachie-Landshoff (DL) with leading twist only, and with higher twist (HT) 
correction (dot-dashed). 
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FIG. 3. Ratio of nuclear EMC ratios for 3 He and 3 H for the PEST wave function (solid), modi- 
fied PEST to reproduce the experimental binding energies (dashed), and the density extrapolation 
model (dot-dashed). 
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FIG. 4. F2/F2 ratio extracted from previous deep- inelastic p and d cross sections using a 
Fermi-smearing model || (solid circles). The shaded band represents a ± one standard deviation 
error for the proposed 3 H and 3 He DIS JLab experiment, with the central values of the band 
corresponding to F^/F^ extracted assuming an EMC effect in deuterium (see text). 
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